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IN THE SPECIFICATION: 

Please amend te specification as follows: 

Please repl ace the originally filed paragraph beginning at page 32> line 7 as follows: 

Figure 5 dem o nstrat e s Figures Sa-Si demonstrate how the ITD may be applied to 
determine a TFE distribution of a sample signal and compares its performance to prior art 
methods of Fourier analysis and wavelet analysis for this task. Figure 5A 5a contains a sample 
signal. Figures 5B=e 5b-5c contain the TFE distribution of the signal obtained via windowed 
Fourier analysis. Note the rectangular grid of TFE information naturally derived fi-om this 
transform and the trade-offs between good temporal but poor frequency locahzation (as shown in 
Figure 5B 5b ), and poor temporal but good frequency localization (as shown in Figure 5€ 5c ). 
Darker regions correspond to those with higher power throughout this Figure. It is impossible to 
simultaneously improve both temporal and fi*equency locahzation of energy. Quantification is 
performed in this Fourier analysis by sequential time windows with a predetermined 
segmentation of the signal that is not determined in any way by the signal changes (e.g., local 
extrema). Figure 5B 5d illustrates the TFE distribution for the same input signal that is obtained 
using the fast wavelet transfor m, sometimes referred to herein as the FWT , Note the dyadic grid 
of TFE information naturally derived fi-om this transform. The FWT has the ability to locaUze 
higher fi-equency information on shorter time scales and lower frequency information on longer 
timescales, but still suffers from predetermined temporal segmentation and inaccuracies in both 
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temporal and frequency localization. These are due in part to the bleeding of signal energy 
across different levels of resolution in the wavelet transform and are attributed in large part to the 
temporal and frequency bins that are predetermined by the choice of basis and not by the signal 
under study. Figure 5E 5e shows the ITD-based TFE distribution for the same signal, illustrating 
the significant advance in time-frequency localization provided by the algorithm of the subject 
invention. Each displayed line segment's start and end points correspond exactly to the start and 
end of an actual wave in an ITD proper rotation component and the line segment shading is 
determined by the amplitude of the corresponding wave (darker indicates greater amplitude). 
Since the ITD automatically decomposes the input signal into a set of proper rotations and a 
monotonic trend, the instantaneous phase angles for each proper rotation component can be 
obtained via the Hilbert transform or the invention's alternative approaches described herein to 
obtain the corresponding instantaneous frequencies. Figure 5F 5f (pagc2of3) shows the 
original signal and the first four HF components obtained by the ITD. Figure 5B 5g shows the 
instantaneous frequency curves (in light gray when instantaneous power is insignificant, i.e., 
below 2x10"^) corresponding to each of the HF proper rotation components (as shown smoothed 
and solid black) when instantaneous power is non-negligible, i.e., above 0.005, to illustrate this 
process. The curves in Figure 5G 5g were obtained using the arcsine-based method for 
instantaneous frequency computation and were smoothed using a 0.1 second median filter to 
enhance visualization. Figure 5H 5h (page 3 of 3) shows the first three proper rotation 
components obtained by application of the ITD to the signal of Figure 5A 5a . Figui ' c 51 Figures 
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5i and 5} 5x show the half-wave amplitudes and (smoothed) instantaneous frequencies for each 
component obtained via the arcsine-based method for instantaneous phase and instantaneous 
frequency detemiination (linestyles are the same as that of the corresponding proper rotation 
components). As expected, the instantaneous frequencies of each component decrease at each 
point in time as the level increases, but the actual frequency values are determined solely by the 
signal and not predetermined in any manner. — 

Please replace the originallv filed paragraph beginning at page 34. line 5 as follows: 

The abiHty of the ITD to decompose online any input signal into a set of one or more 
component signals which have the proper rotation property makes it an ideal first step of a two- 
step process resulting in a powerfiil and unique new method for nonlinear signal filtering. In 
particular, as each monotonic segment is computed within one or more decomposition levels, 
features of the segment can be analyzed and quantified. For example, the amplitude and duration 
of the segment may be simply computed from the extrema, i.e., the starting and ending points of 
the segment. One skilled in the art will appreciate that more detailed features such as the time- 
weighted feature density of the segment , sec for example hit'l. Publication No. WO 01/75660 
Al, o fNikitin c t al, and associated properties such as its average or median value, its variance, 
skewness, distance from a template or reference density, etc., may also be quantified for each 
segment. Individual segments may then be classified according to the values of theses quantified 
features. A filtered signal can then be constructed by adding together (i.e., superimposing the 
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segments while preserving the time intervals on which each occurred) only those segments that 
satisfy certain constraints applied to their quantified features. For example, one may construct a 
signal using only those ITD-created monotonic segments with durations that, using the data 
sampling rate, correspond to wavespeeds between Hz and 4,ax Hz, and with absolute 
amplitudes of at least A^j^, e.g., in excess of the 75* percentile of amplitudes for all segments 
with the required duration. In this example, the output can be essentially interpreted as extracting 
the waves in the input signal that are in the f^^ - f^^^x frequency band and which have larger 
amplitudes than A„i„. Such uses would bear a strong resemblance to the commonly used Fourier- 
based band-pass filtering, but would avoid the phase-shifting and waveform distortion 
drawbacks, while adding the amplitude differentiation capability. In Figure 6 Figures 6a-6g . this 
concept is illustrated when applied to the same input signal used in Figure 5A 5a (also shown 
again in Figure 6A 6a ). After the TFE distribution is computed, the absolute amplitude and 
duration/wavespeed of each individual wave are obtained. Figure 6D-C shows Figures 6b-6c 
show the one-dimensional (marginal) densities of these two features along with vertical 
discriminator in solid lines that were selected for purposes of illustrating the capabilities of this 
method. Figure 6B 6d shows the (absolute amplitude, wavespeed) ordered pairs for each 
segment. Discriminators of pairs with wavespeed less than 20 Hz, between 20 Hz and 50 Hz, 
and greater than 50 Hz are applied, along with a discriminator of amplitude exceeding 0.01 mV. 
Feature vector pairs (i.e., points in Figure ^ 6d ) falling into certain specific regions of the two- 
dimensional feature range are classified according to the region in which they belong. Figure 6E 
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6e displays three filtered components constructed using only those segments with points in 
specific regions (corresponding to three of the regions in Figure 6B 6d ), Component one, Cj, is 
the output of a filter that contains the signal reconstructed using only those waves that have 
absolute amplitude exceeding O.OlmV and wavespeed between 20Hz and 50 Hz. Component 
two, C2, is the output of a fiher that contains the signal reconstructed using only those waves that 
have absolute ampHtude less than O.OlmV and wavespeed between 20Hz and 50 Hz. Component 
three, C3, is the output of a filter that contains the signal reconstructed using only those waves 
that have absolute amplitude exceeding O.OlmV and wavespeed below 20Hz. The output 
obtained by combining some of these components (i.e., C1+C2+C3 and Cj+Cj) are shown in 
Figure €F6f, along with the output of a traditional, linear, [20, 50] Hz. band-pass digital filter 
provided for comparison. The extrema-preserving properties of this new type of filter are evident 
in this Figure 6F 6f. Figure 66 6g again compares an output available by using the ITD filter in 
comparison to the same linear band-pass digital filter, but does so on a sUghtly larger timescale 
(10 seconds). The residual signals obtained by subtracting the filtered outputs firom the original 
signal are also shown, to illustrate the ability of the ITD based signal to more completely 
decompose the raw input signal into what may be interpreted as abnormal and normal 
components present in the raw signal. - 
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Please replace the originally filed paragraph beginning at page 36. line 23 as follows: 

Another example that illustrates the power of the invention for signal filtering and 
analysis is provided in Figui ' c 7 Figures 7a-7h . In tills Figui - c these Figures, we demonstrate how 
the ITD-based filtering method may be applied to differentiate between two types of waves that 
have significantly overlapping spectral characteristics and to decompose a raw signal that is a 
combination of these two types of waves into components that preserve temporal location of 
extrema and critical points in the raw signal. Overlap between the power spectral densities of an 
underlying signal and superimposed noise is one of the more difficult problems in signal analysis 
and detection, with numerous applications. Figure ?A 7a illustrates the mix of two altemating 
signals, one a 9 Hz cosine wave with an increasing linear trend and the other a 9 Hz sawtooth 
wave. The signals were added to create a raw test signal for use in this example. Figure ?B 7b 
shows the power spectral density estimates for each of the two signals, clearly illustrating the 
significant overlap in the band around 9 Hz. Figure 9G 7c shows the proper rotation components 
and monotonic trend obtained via appUcation of ITD to the original combined signal. Figure 7D- 
B Figures 7d-7e illustrate the result of applying single wave analysis to examine various features 
of individual waves obtained by the ITD decomposition. Figure 7d shows the ratio of 
individual wave peak to mean signal value, while Figure 7e shows a plot of wave skewness 
plotted versus wave kurtosis. A sample discriminator in dashed linestyle is also shown to 
illustrate the result of feature classification as applied to these quantified waveform features. 
Figure ?F 7f shows the output of the ITD-based filter that retains only those waves with peak-to- 
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mean value ratio exceeding 1.65. The corresponding residual signal, i.e., the original signal 
minus the filter output signal, is also shown and illustrates that peak-to-mean value ratio 
discrimination of ITD-produced proper rotation waves allow reconstruction of a filter output that 
nearly perfectly decomposes the original signal into its constituent components while retaining 
precise temporal information regarding, e.g., extrema location. Figure 7g shows the output 
of the ITD-based filter that retains only those waves with skewness-kurtosis pairs above the 
dashed discriminator line, along with the corresponding residual signal. This example illustrates 
how this new filter method may be easily configured to decompose an input signal into 
components with more complex criteria, such as "cutting off the tops of all riding sinusoidal 
waveforms while retaining every other signal component" as is done here. Figure 9H 7h shows 
the same signals as in Figure 7G but zooming in to the time interval between t=9 and t=13 to see 
details of the filter output.-- 

Please replace the originally filed paragraph beginning at page 39. line 5 as follows: 

One skilled in the art will also recognize that, just as Fourier-based filtering and other 
decomposition methods can be successfiiUy applied to higher dimensional analysis, so, too, can 
the ITD system. This extends the above-mentioned appUcation of the ITD to decompose 
functions of several variables. For example, in image processing applications, one may be 
interested in separating an image into a high firequency component, which contains edge 
transitions between various objects in the image for example, and a low fi-equency component, 
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which define background colors in the image for example. Two-dimensional wavelet and/or 
Fourier analysis are popular for image processing applications. However, the drawbacks 
mentioned above for these methods still exist in higher dimensional analysis. The benefits of the 
ITD system, such as its ability to preserve extrema and precisely localize time-fi-equency-energy 
information are preserved in higher dimensional analysis and thereby offer improvements over 
the most popular prior art methods. Figure 8 illu&Ualcs Figures 8a-8c illustrate the application of 
ITD to decompose a two-dimensional surface (shown in Figure ftA 8a) into two component 
surfaces, namely a high frequency surface (shown in Figure SB 8b) and a low frequency surface 
(shown in Figure 8c). The decomposition in this example is obtained by decomposing each 
of the cross-sectional signals of the surface along the grid lines obtained by holding the first 
independent variable constant, repeating the process holding the second variable constant, and 
averaging the two results. - 

After the originally fil ed description of Figure 8 on page 15. please insert the following : 

Figure 9a is a flow chart representation of the ITD method of the present invention. 

Figure 9b is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9B of Figure 9a. 

Figure 9C is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9C of Figure 9a. 

Figure 9D is a detailed flow chart representation of a portion of the ITD method of the 
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present invention as indicated at 9D of Figure 9a. 

Figure 9E is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9E of Figure 9a. 

Figure 9F is a detailed flow chart representation of a portion of the ITD method of the 
present invention as indicated at 9F of Figure 9a. ~ 
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